A Planetary Companion to the Intermediate-Mass Giant 

HD 100655 



Masashi Omiya, 1 ' 2 Inwoo Han, 1 Hideyuki Izumiura, 3,4 Byeong-Cheol Lee, 1 Bun'ei Sato, 2 
Kang-Min Kim, 1 Tae Seog Yoon, 5 Eiji Kambe, 3 Michitoshi Yoshida, 6 Seiji Masuda, 7 Eri 

Toyota, 8 Seitaro Urakawa, 9 and Masahide Takada-Hidai, 10 
1 Korea Astronomy and Space Science Institute, 61-1 Whaam-dong, Youseong-gu, Daejeon 305-348, 

South Korea 

2 Department of Earth and Planetary Sciences, Tokyo Institute of Technology, 2-12-1 Ookayama, 

Meguro-ku, Tokyo 152-8551, Japan 
omiya. m@geo. titech. ac.jp 
3 Okayama Astrophysical Observatory, National Astronomical Observatory of Japan, Asakuchi, 

Okayama 719-0232, Japan 
4 Department of Astronomical Science, The Graduate University for Advanced Studies, Shonan 

Village, Hayama, Kanagawa 240-0193, Japan 
5 Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Daegu 

102-101, South Korea 

6 Hiroshima Astrophysical Science Center, Hiroshima University, Higashi-Hiroshima, Hiroshima 

139-8526, Japan 

7 Tokushima Science Museum, Asutamu Land Tokushima, Itano-gun, Tokushima 119-0111, Japan 
8 Kobe Science Museum, 1-1-6 Minatojimanakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan 
9 Bisei Spaceguard Center, Japan Spaceguard Association, 1116-3 Okura, Bisei-cho, Ibara, Okayama 

114-1411, Japan 

10 Liberal Arts Education Center, Tokai University, 4-1-1 Kitakaname, Hiratsuka, Kanagawa 

259-1292, Japan 

(Received ; accepted ) 
Abstract 

A precise radial velocity survey conducted by a Korean— Japanese planet search 
program revealed a planetary companion around the intermediate-mass clump giant 
HD 100655. The radial velocity of the star exhibits a periodic Keplerian variation 
with a period, semi-amplitude and eccentricity of 157.57 d, 35.2 m s _1 and 0.085, 
respectively. Adopting an estimated stellar mass of 2.4 M Q , we confirmed the presence 
of a planetary companion with a semi- major axis of 0.76 AU and a minimum mass 
of 1.7 Mj. The planet is the lowest-mass planet yet discovered around clump giants 
with masses greater than 1.9 M Q . 
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1. Introduction 

Over 550 exoplanets have been discovered to date. Many of the planets orbit solar- 
mass (0.7—1.5 M Q ) stars, and they have revealed properties that are now used to constrain 
planet-formation models (e.g., Ida & Lin 2004; Butler et al. 2006; Udry & Santos 2007). In 
contrast, only about 60 and 25 planets have been detected around evolved G-K (sub)giants 
(1.5—5 M Q ) and K-M dwarfs (<0.7 M ), respectively (e.g., Sato et al. 2008; Johnson et al. 
2011; Johnson et al. 2007a). Accordingly, the properties of the planetary systems orbiting such 
stars are less clarified yet than those for solar-mass stars. Planetary formation depends on 
the properties of protoplanetary disks, which should be affected by properties of the host star, 
such as stellar metallicity, radiation output, and disk diffusion times (e.g., Kornet et al. 2006; 
Kennedy & Kenyon 2008). Observational features of planetary systems over a wide range of 
host star masses need to be clarified by current and future surveys of various masses stars in 
order to understand planetary formation in general. 

More than 20 years ago, initial theoretical ideas of planetary formation for systems over a 
wide range of stellar masses were presented in terms of planet formation in protoplanetary disks 
with different properties (Nakano 1988a; Nakano 1988b). In the last two decades, improvements 
in planet formation modeling have made it possible to compare theoretical models directly with 
observed properties of planetary systems around stars with various masses (e.g., Ida & Lin 
2005; Burkert & Ida 2007). For example, Kennedy & Kenyon (2008) predicted that the peak 
occurrence rate of giant planets occurs for stars with masses of around 3 M , based on a core 
accretion scenario which includes the movements of snow lines under the evolution of central 
stars. Moreover, Currie (2009) suggested that "the planet desert", i.e., a dearth of planets 
with semi-major axes of <0.6 AU orbiting >1.5 M & stars, may be reproduced by the effects of 
Type-II migration, considering the dependence of diffusion time of the protoplanetary disk on 
stellar mass. Clarifying the relationship between stellar mass and planetary system will provide 
valuable insights into planet formation models. 

For intermediate-mass stars on the main sequence, precise Doppler surveys are difficult 
because of their large intrinsic radial velocity variations and smooth spectra with few absorption 
lines, caused by high surface activity, high surface temperature and/or high rotational velocity 
(Lagrange et al. 2009). In contrast, evolved intermediate-mass (sub)giant stars are suitable 
targets for precise Doppler surveys because these stars have low surface activity and their 
spectra exhibit many sharp absorption lines. Thus, to date, spectroscopy-based planet searches 
targeting intermediate-mass stars have been carried out through precise Doppler surveys of 
evolved stars. Although the number of substellar companions found orbiting such stars is still 
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insufficient, some characteristic planetary system properties across a wide range of host star 
masses have begun to emerge. For example, the masses of planets and their host stars show 
correlation: more massive substellar companions tend to exist around more massive stars (e.g., 
Lovis & Mayor 2007). This correlation suggests that the mass range of the brown dwarf desert 
depends on host-star's mass, and that planets may be deficient around 2.4—4 M & stars (Omiya 
et al. 2009). Also, the planet occurrence rate depends on host-star's mass: the giant planet 
frequency for higher-mass giant stars is higher than that for lower-mass stars (Lovis & Mayor 
2007; Johnson et al. 2007a). The fraction of giant planets increases with increasing stellar mass 
up to 2 M & (Johnson et al. 2010a). Moreover, the orbital semi-major axes of planetary systems 
also seem to be correlated to host-star's properties. Semi-major axes of most planets orbiting 
intermediate-mass (sub)giant stars are larger than 0.6 AU 1 , while those orbiting solar-type stars 
are larger than 0.02 AU (Johnson et al. 2007b; Sato et al. 2008; Wright et al. 2009, Bowler 
et al. 2010). Even considering the effect of engulfment of inner-orbit planets by host stars, 
which have experienced rapid expansion in the red giant branch (RGB) phase (Sato et al. 2008, 
Villaver & Livio 2009), the observed properties of substellar systems orbiting intermediate- mass 
(sub)giant stars seem to be different from those orbiting solar-type stars (see also Bowler et al. 
2010). 

In 2005, we started a Doppler spectroscopy-based survey of evolved GK-type giants in 
a framework of a Korean— Japanese planet search program (Omiya et al. 2009). The survey 
program is an extension to the ongoing Okayama Astrophysical Observatory (OAO) planet 
search program (Sato et al. 2005), and aims to clarify the properties of their associated planetary 
systems in collaboration with an East-Asian Planet Search Network (EAPS-Net; Izumiura 
2005). About 190 sample stars of the survey were selected from the Hipparcos catalog based 
on the same criteria as those for OAO planet search program, except visual magnitude (6.2 
< V < 6.5). The radial velocity variability of each sample star is monitored using either the 
1.8-m telescope at Bohyunsan Optical Astronomy Observatory (BOAO, Korea) or the 1.88-m 
telescope at OAO (Japan). If a sample star exhibits large variations in radial velocity, follow-up 
observations of the star are performed using both telescopes. 

In this paper, we report the discovery of a planetary companion orbiting the 
intermediate-mass giant HD 100655. This is the first planet discovered by this Korean- Japanese 
planet search program. In section 2, we describe our observations and radial velocity measure- 
ments from BOAO and OAO data. The properties of the host star and the radial velocity 
variability are reported in sections 3 and 4, respectively. We discuss possible causes of the 
radial velocity variation in section 5. In section 6, we consider the implications of this discovery 
for the current picture of planetary companions around intermediate-mass giant stars. 

1 A planet with a semi-major axis of 0.081 AU was found orbiting an intermediate-mass subgiant star HD 
102956 with a mass of 1.68 M© (Johnson et al. 2010b). 
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2. Observations and Analyses 

2.1. BOES Observations and Analysis 

Radial velocity observations at BOAO were carried out with the 1.8-m telescope and 
the BOAO Echelle Spectrograph (BOES; Kim et al. 2007), a fiber-fed high resolution echelle 
spectrograph. We placed an iodine (I2) cell in the optical path in front of the fiber entrance of 
the spectrograph (Kim et al. 2002) for precise wavelength calibration and used a 200-/iin fiber, 
obtaining a wavelength resolution R = A/AA ~ 51,000. The spectra covered a wavelength 
region from 3500 A to 10,500 A. Echelle data reduction was performed using the IRAF 2 
software package in the standard manner. We used a wavelength region of 5000—5900 A which 
is covered by many I2 absorption lines, for radial velocity measurements. We also made use of 
Ca II H line at around 3970 A as chromospheric activity diagnostics. Radial velocity analysis 
was performed using the spectral modeling technique described in Sato et al. (2002), which was 
based on the method of Butler et al. (1996) and was adapted to BOES data analysis (Omiya 
et al. 2009). We employed the extraction method described in Sato et al. (2002) to prepare a 
stellar template spectrum from stellar spectra taken through the I2 cell (star+12 spectra). The 
technique allowed us to achieve a long-term Doppler precision of 14 m s" 1 over 4.5 years. 

2.2. HIDES Observations and Analysis 

Radial velocity observations at OAO were carried out with the 1.88-m telescope and 
High Dispersion Echelle Spectrograph (HIDES; Izumiura 1999) attached to the coude focus of 
the telescope. We used an I 2 cell placed in the optical path in front of the slit of the spectrograph 
(Kambe et al. 2002) as a precise wavelength calibrator. We always set the slit width to 200 
/im (0.76"), providing a spectral resolution of 63,000. Until November 2007, we had taken 
star+I 2 spectra with a wavelength region of 5000— 6200-A. Since the HIDES CCD system was 
upgraded to a three-CCD mosaic in December 2007, we have obtained spectra from 3750 A 
to 7550 A. The wavelength region of 5000—5900 A of the star+I 2 spectra are used for radial 
velocity measurements. The full range of stellar spectra taken without the I2 cell are used for 
abundance analysis. Echelle data reduction was performed using the IRAF software package 
in the standard manner. Stellar radial velocities were derived from the star+I 2 spectra using 
the spectral modeling techniques detailed in Sato et al. (2002), giving a Doppler precision of 
less than 8 m s" 1 over 4.5 years. 

2 IRAF are distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc. under cooperative agreement with the National 
Science Foundation, USA. 



4 



3. Stellar Parameters of HD 100655 



HD 100655 (HR 4459, HIP 56508, BD+21 2331) is 122.3 ± 7.5 pc from the Sun according 
to the Hipparcos parallax of n = 8.18 ± 0.50 mas (van Leeuwen 2007). The star is classified 
as a G9III giant star with V = 6.45 and B - V = 1.010 ± 0.015 (ESA 1997). We corrected 
the observed color index by an extinction value of EiB — V) = 0.0163 ± 0.0016. The value 
was calculated from the galactic extinction of E(B — V)$ = 0.0273 ± 0.0015 to the direction of 
the star obtained from the Schlegel, Finkbeiner and Davis (1998) dust maps using the relation 
E(B — V) — E(B — V)s[l— exp(— |JDsinb|/125)], where D and b are the distance from sun and 
the galactic latitude, respectively. We derived an effective temperature of the star of T e g = 4861 
± 110 K using the (B — V) — T cS calibration of Alonso et al. (1999, 2001). A luminosity of L 
= 43 ± 5 L Q was obtained from the absolute magnitude My = 0.96 ± 0.13 and the bolometric 
correction B.C. = —0.31 ± 0.04 based on the calibration of Alonso et al. (1999). A stellar mass 
of M = 2.4±8;1 M© was estimated by interpolating the evolutionary tracks of Girardi et al. 
(2000) with the estimated T eff and L (see figure 1). We determined the surface gravity to be log 
g = 2.89 ± 0.10 and the stellar radius R = 9.3t\'_i Rq from M, L, and T e ff- The microturbulent 
velocity V t = 1.36 ± 0.03 km s -1 and the [Fe/H] of 0.15 ± 0.12 were derived from abundance 
analysis of a model atmosphere (Kurucz 1993) using the equivalent widths of Fe I and Fe II 
lines measured from an l2-free spectrum of HD 100655. We adopted gf- values of Fe I and Fe 
II lines from Takeda et al. (2005). de Medeiros & Mayor (1999) found the stellar rotational 
velocity, vsini s , to be 1.6 ± 1.0 km s -1 . This value is comparable to the rotational velocities of 
typical late G-type giants. The stellar parameters are summarized in table 1. 

4. Orbital Solution 

A large radial velocity variation in the star HD 100655 was found in the early BOAO 
survey and we made intensive follow-up observations of the star at BOAO and OAO For 4.5 
years from the beginning of the survey, we collected 13 BOAO data points having a typical 
signal-to-noise ratio (S/N) of 170 pixel -1 with an exposure time of 900—1200 s, and 32 OAO 
data points having a typical S/N of 120 pixel -1 with an exposure time of 1200—1800 s. The 
observed radial velocities of HD 100655 are shown in figure 2 and listed in table 2, together with 
the observation dates (JD) and estimated uncertainties. A dominant peak in the Lomb-Scargle 
periodogram (Scargle 1982) of the radial velocity variation exists at a period of 157.78 d (a 
frequency of 0.006338 c d -1 ) (see figure 3). To check the significance of this periodicity, we 
estimated a False Alarm Probability (FAP) using the bootstrap randomization method. We 
produced 10 5 fake data sets by randomly mixing the observed radial velocities with a fixed 
observation date, and applied the Lomb-Scargle periodogram analysis to them. Only one fake 
data set showed a periodogram power higher than the observed one. Thus, the FAP of the 
period is 10 -5 . A best-fit Keplerian orbit derived from both the BOAO and OAO velocity data 
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by a least-squares fit has a period P = 157.57 d, a velocity semi-amplitude K\ = 35.2 m s , 
and an eccentricity e = 0.085. The best-fit curve is shown in figure 2 as a solid line overlaid 
on the observed velocities. We applied an offset of ARV = —28.1 m s" 1 to the BOAO velocity 
data, estimated concurrently with the orbital fit to a Keplerian model. The offset was required 
because of difference of velocity zero points between BOES and HIDES data originated from 
using different stellar templates for each data. The rms of the residuals to the best-fit are 14.9 
m s" 1 for BOAO data, 9.2 m s _1 for OAO data, and 11.2 m s _1 for combined data sets. In the 
residuals we could not find any significant periodic variation due to additional companions. The 
best-fit orbital parameters and their uncertainties are listed in table 3. The uncertainties were 
estimated using a bootstrap Monte Carlo approach by creating 1000 fake data sets. Adopting 
a stellar mass M = 2.4+o 2 4 M & for HD 100655, we obtained a semi-major axis a — 0.76l(xo4 
AU and a minimum mass A^sim^ = 1.7lo]2 f° r the planetary companion. 

5. Cause of the Radial Velocity Variation 

To examine causes of the apparent radial velocity variation other than orbital motion, we 
checked the Ca II H line and the Hipparcos photometric variation, and performed spectral-line 
shape analyses using a technique described in Sato et al. (2007) as follows. In the analyses, we 
investigated the cause of the velocity difference between spectra observed at top and bottom 
velocity phase. 

Figure 4 shows the spectrum around the Ca II H line of HD 100655. We note a lack of 
significant emission in the Ca II H line core of HD 100655, which suggests chromospheric in- 
activity for the star. Moreover, Hipparcos photometry demonstrates the photometric stability 
of HD 100655 down to a ~ 0.008 mag. based on the 55 observations for the star over a period 
of 1000 d. Figure 5 displays a periodogram of the Hipparcos photometry. We note a weak 
peak around the period of the radial velocity variation. To check the significance of the peak, 
we estimated FAP using the bootstrap method as well as the method described in section 4. 
We produced 10 5 fake data sets, and applied the Lomb-Scargle periodogram analysis to them. 
A total of 7425 fake datasets showed a peak around the period of the radial velocity variation 
higher than the peak on the observed data set, which means FAP of the peak is about 7.4 
%. Thus the peak is not considered to be significant. Although we have not completely dis- 
proved the possibility that the radial velocity variation is due to rotational modulation, these 
photometric results suggest that the main cause of the observed radial velocity variation is not 
rotational modulation of stellar spots. 

For spectral-line shape analysis, we extracted two high-resolution stellar templates from 
star(HD 100655)+i2 spectra obtained at OAO, using the method described in Sato et al. (2002). 
One template was constructed from four spectra with observed radial velocities of the peak 
phase (~32 m s _1 ), and the other from four spectra of the valley phase (—44 m s" 1 to —34 
m s" 1 ). Cross-correlation profiles of the two templates were provided for 75 spectral segments 
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(4- A to 5-A width each) that did not include severely blended lines or broad lines. We obtained 
a bisector for the cross-correlation profile of each segment and calculated three quantities from 
velocities at three flux levels (25%, 50%, and 75%) of the bisector profile. One quantity is 
the bisector velocity span (BVS), which is the velocity difference between two flux levels with 
25% and 75% of the bisector. Another is the bisector velocity curvature (BVC), which is the 
difference between two velocity spans in the upper half (between two flux levels with 50% and 
75% of the bisector) and the lower half (25% and 50%). The other is the bisector velocity 
displacement (BVD), which is the average of the velocities at the three flux levels (25%, 50%, 
and 75%). These bisector quantities for HD 100655 are shown in figure 6. The average values 
of BVS and BVC are 7.8 ± 8.1 m s _1 and 2.4 ± 3.9 m s _1 , respectively. The BVS values may 
be increased due to rotational stellar spots, which may invoke photometric variation. However, 
since the average value of the BVS is one ninth of the velocity differences (~70 m s" 1 ) between 
the two templates, we consider both BVS and BVC value to be essentially zero, meaning that 
the cross-correlation profiles are symmetric. Moreover, the average value of the BVD (—70.1 
± 17.9 m s _1 ) is consistent with the velocity difference between the two templates. Thus, 
the cause of the velocity difference is considered to be a parallel shift of spectral lines, not 
variations in spectral line shapes. Hence, the observed radial velocity variation of HD 100655 is 
best explained by the orbital motion of a planetary companion, not by intrinsic activity, such 
as rotational modulation and pulsation. 

6. Discussion 

We detected a planetary companion orbiting the clump giant star HD 100655 based 
on the precise Doppler spectroscopy survey conducted by the Korean— Japanese planet search 
program. The radial velocity variation of the star discovered during early observation at BOAO 
indicated the existence of a possible planetary companion, and the orbital parameters of the 
companion were determined by follow-up observations at BOAO and OAO. Adopting a mass 
of 2.4 M Q for HD 100655, we found that the planetary companion has a minimum mass of 1.7 
Mj and a semi- major axis of 0.76 AU. This is the lowest-mass planet among those discovered 
around giant stars with masses larger than 1.9 M . Fourteen planetary companions and six 
brown dwarf-mass companions have been detected so far around such giants by ongoing precise 
Doppler surveys, and the discoveries bring out some characteristic properties of the planetary 
systems. 

Figure 7 plots mass of substellar companions with semi-major axes less than 3 AU 
against host-star's mass (updated version of figure 5 of Omiya et al. 2009). This figure includes 
intermediate-mass (1.5 M Q < M < 4 M & ) giants and subgiants (filled circles), intermediate- 
mass dwarfs (open circles), solar-mass stars (M < 1.5 M & , open triangles), and HD 100655 
(star). Solid and dot-dashed lines indicate the lower-mass limits of companions detectable by 
current Doppler surveys for semi-major axes of 0.6 AU and 3 AU, respectively. The limits 
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correspond to companion masses that give rise to semi-amplitudes of radial velocity variations 
of their host stars as large as three times the typical radial velocity jitters, which are 5 m s" 1 for 
subgiants (1.5—1.9 M & ) and 20 m s" 1 for clump giants (1.9—4 M & ) (Johnson et al. 2010c, Sato 
et al. 2005). Two unpopulated regions of substellar companions orbiting intermediate-mass 
subgiants and giants appear in regions (a) and (b) in figure 7 (Omiya et al. 2009). The planet 
orbiting HD 100655 is located below the detection limit and the region (b) in figure 7 because 
of its small root mean-square scatter of the residual radial velocities (~11 m s" 1 ), that is, its 
small radial velocity jitter. The existence of this planet suggests a possibility that low- mass 
giant planets can form around ~2.4 M & stars, while this planet could possibly have a small 
orbital inclination, and thus a high actual mass. Therefore, a paucity of low-mass companions 
orbiting massive intermediate-mass giants, roughly indicated by the region (b), might partly be 
caused by an observational bias due to the high detection limit. In this respect, observational 
surveys more sensitive to lower-mass substellar companions are necessary. 

The mass distribution of substellar companions orbiting 1.5—3 M stars may also depend 
on the semi-major axes of the companions. Figure 8 is a plot of semi-major axis of substellar 
companion versus host star mass. Crosses, circles and filled circles indicate brown dwarf-mass 
companions (13—30 Mj), " superplanets" (6—13 Mj), and normal giant planets (1—6 Mj), 
respectively. Solid, dot-dashed and dotted lines indicate the typical farthest orbital distances 
of companions detectable by current Doppler surveys for companion masses of 3, 4 and 5 Mj, 
respectively. The distances correspond to orbital semi-major axes that the companions induce 
radial velocity variations of their host stars with semi-amplitudes as large as three times the 
typical radial velocity jitter, which is 20 m s" 1 for clump giants (1.9—3 M Q ) (Sato et al. 2005). 
In figure 8, some interesting properties of substellar companions are suggested in three stellar 
mass ranges. Almost all the planets orbiting 1.5—1.9 M stars are normal giant planets, and 
are located on orbits with semi-major axes of >1 AU. Many planets orbiting 1.9—2.5 M & stars 
seems to be classified in two groups 3 : normal giant planets at inner orbits (0.6—1.3 AU) and 
superplanets at outer orbits (1.9—3 AU). HD 100655 b is included in the group of the normal 
giant planets. All planet-mass companions orbiting 2.5—3 M Q stars reside at semi-major axes 
larger than 1.9 AU, while all brown dwarf-mass companions are orbiting at semi-major axes less 
than 1.9 AU. Although the number of known substellar companions discovered around stars 
with >2.5 M Q is still small, the distribution of substellar companions around 1.9—2.5 M & stars 
may differ from those around 1.5—1.9 M & and 2.5—3 M stars. 

To reproduce the distribution of giant planets around 1.9—2.5 M Q giant stars, two 
scenarios can be suggested. One is the planet engulfment scenario caused by stellar evolution 
of primary stars. Most of the host stars are clump giants that should have experienced the 
RGB phase, which triggers rapid stellar expansion. Villaver & Livio (2009) suggested that the 

3 We note that <3 Mj (<5 Mj) planets orbiting such stars at 1 AU (3 AU) are below the lower-mass limits 
for typical detectable planets. 
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primary stars can preferentially capture more massive planetary companions by tidal interaction 
in the RGB phase. Thus, the superplanets with semi-major axes of <1.9 AU might have 
been preferentially engulfed by their primary stars even if they had existed, leaving normal 
giant planets. However, according to Kunitomo et al. (2011), the critical semi-major axis, 
within which a primary star can engulf planetary companions, decreases from ~ 1.5 (0.4) AU for 
1.7(2.0) M Q stars to ~0.2 AU for 2.1 M stars and thus all the planets with semi-major axis 
larger than 0.6 AU around 2.0—2.5 M stars can survive the RGB phase regardless of their 
masses. Therefore, the observational properties would not be quantitatively explained by only 
this mechanism. 

The other scenario is that the distribution of the planets is primordially originated from 
planet migration in protoplanetary disk. Dependence of Type-II migration rate on planet mass 
may separate locations of low-mass giant planets and superplanets. For example, based on 
the equation (1) of Currie (2009), 2 Mj and 8 Mj planets that formed in circular orbits with 
a semi-major axis of 3.8 AU around 2 M Q stars can migrate to inner orbits with semi-major 
axes of ~0.7 AU and ~2.8 AU, respectively, assuming a disk dissipation time of 1 Myr. Thus, 
the observed orbital distribution of planets around 1.9—2.5 M Q stars may be explained by this 
mechanism, and if this is the case many undetected lower-mass planets should be expected at 
distances of 1—3 AU because giant planets can form at any distance beyond the snow line. It 
should be noticed, however, the observed semi-major axis distributions of planetary systems 
around 1.5—1.9 M Q and 2.5—3 M & stars might not be explained by only the effect of the 
migration. 

Additionally, Type-II migration may not be only the mechanism that can locate gi- 
ant planets around intermediate-mass stars. The magnetorotational instability-dead zone in 
the protoplanetary disks may encourage formations of giant planets only at ~1 AU around 
intermediate-mass stars (Kretke et al. 2009). In this case, a drop-off of giant planets at >~1 
AU might exist around such stars. However, the distribution of <3 Mj planets at larger than 
1 AU around 1.9—3 M & stars has not been clarified yet due to the detection limits of current 
planet searches. 

Thus, in order to examine roles of these mechanisms on planet formation and evolution 
around intermediate-mass stars, it is required to evaluate the semi-major axes distribution 
by further Doppler surveys of intermediate- mass stars with masses of >1.9 M sensitive to 
lower-mass planets. 
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Fig. 2. Upper panel: radial velocities of HD 100655 observed at BOAO (filled circles) and OAO (open 
circles). The solid line represents the Keplerian orbital curve. Lower panel: Residuals to the best Keplerian 
fit. 
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Fig. 3. The Lomb-Scargle periodgram of the radial velocity variation of HD 100655. A dominant peak 
appears at a period of 157.78 d (a frequency of 0.006338 c d _1 ) with a False Alarm Probability (FAP) of 
10~ 5 . 
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Fig. 4. The spectrum around the HD 100655 Ca II H line. The line core does not seem to exhibit high 
chromosphcric activity. 
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Fig. 5. Pcriodgram of Hipparcos photometric variation in HD 100655. Although a weak peak around 
the period of the stellar radial velocity variation does appear (arrowed line), a FAP of the peak is about 
7.4 %; thus, the peak is not significant one. 
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Fig. 6. Bisector quantities obtained from calculations of cross-correlation functions of two distinct stellar 
templates. The templates are constructed from star+I 2 spectra with radial velocities of the peak and valley 
phase. Values of bisector velocity span (BVS, circles), bisector velocity curvature (BVC, triangles) , and 
bisector velocity displacement (BVD, squares) are shown with offsets of 400 m s _1 , m s _1 and —500 m 
s _1 , respectively, and their offsets are represented by the dotted- lines. 




Fig. 7. Planetary mass and stellar mass of the planetary systems. Solid dots, circles, and triangles 
represent planetary systems orbiting intermediate-mass (sub)giants, intermediate-mass dwarfs and solar 
mass stars, respectively. A star indicates the planetary system orbiting HD 100655. Dot-dashed and solid 
lines mark the detection limits of companions orbiting stars of any mass at 0.6 and 3 AU. Two of the 
unpopulated regions shown by Omiya ct al. (2009) are indicated by (a) and (b). The planetary system of 
HD100655 appears below the detection limit, because of its small orbital semi-major axis and low radial 
velocity jitter. 
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Fig. 8. Orbital semi-major axis versus stellar mass of planetary systems. Dots, circles, and crosses indi- 
cate the locations of giant planets (1—6 Mj), supcrplanets (6—13 Mj), and brown dwarfs (13—30 Mj), 
respectively. Solid, dot-dashed and dotted lines indicate the typical largest semi-major axes of companions 
detectable by current Doppler surveys with masses of 3, 4 and 5 Mj, respectively. Many planets around 
1.9—2.5 Mq giant stars seem to belong to normal giant planets orbiting at semi-major axes of 0.6—1.3 
AU, or superplanets orbiting at semi-major axes of 1.9—3 AU. 



Table 1. Stellar parameters of HD 100655 

Parameter Value 
Spectral Type 
V 

B-V 
7r (mas) 

My 

B.C. 

Tcff (K) 
L (L ) 

M (Mq) 

R(R e ) 
log 9 

V t (km s- 1 ) 
[Fe/H] 

wsim' s (km s _1 ) 



G9III 
6.45 

1.010 ± 0.015 
8.18 ± 0.50 
0.96 ± 0.13 
-0.31 ± 0.04 
4861 ± 110 
43 ± 5 



2.89 ± 0.10 
1.36 ± 0.03 
0.15 ± 0.12 
1.6 ± 1.0 
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Table 2. Radial velocities of HD 100655 



JD 


Radial Velocity 


Uncertainties 




-2450000 


(m s _1 ) 


(m s^ 1 ) 


Observatory 


3428.2166 


3.9 


8.4 


BOAO 


3809.1443 


77.2 


9.2 


BOAO 


3889.0614 


-4.5 


8.7 


BOAO 


4024.3429 


-12.0 


9.6 


OAO 


4047.3387 


-37.4 


5.6 


OAO 


4075.2415 


2.2 


5.3 


OAO 


4081.3382 


58.1 


10.4 


BOAO 


4094.3171 


8.5 


4.9 


OAO 


4123.2074 


67.9 


8.8 


BOAO 


4131.2972 


30.4 


5.1 


OAO 


4151.2206 


7.7 


4.4 


OAO 


4176.2617 


-8.8 


7.2 


OAO 


4214.1201 


-6.2 


7.1 


OAO 


4224.0873 


33.2 


12.1 


BOAO 


4243.0406 


23.1 


7.5 


OAO 


4262.0463 


16.2 


4.5 


OAO 


4452.3084 


22.6 


9.1 


OAO 


4471.3229 


30.0 


8.1 


BOAO 


4491.2702 


-21.4 


4.3 


OAO 


4505.7574 


-24.3 


8.8 


BOAO 


4527.2340 


-27.8 


4.0 


OAO 


4561.1034 


14.6 


6.4 


OAO 


4565.0837 


18.2 


4.4 


OAO 


4594.0472 


43.3 


5.3 


OAO 


4634.9712 


-6.7 


3.4 


OAO 


4790.3591 


-20.6 


6.3 


OAO 


4816.3500 


-39.9 


4.4 


OAO 


4833.3367 


14.4 


9.0 


BOAO 


4863.2825 


2.9 


7.2 


OAO 


4863.3410 


0.6 


7.4 


OAO 


4881.1750 


20.7 


5.4 


OAO 


4881.3355 


17.3 


4.9 


OAO 


4913.2182 


31.7 


4.0 


OAO 


4927.1213 


34.7 


4.2 


OAO 


4930.0152 


45.3 


10.0 


BOAO 


4931.0597 


31.7 


7.7 


BOAO 


4934.2081 


32.0 


8.0 


OAO 


4937.0928 


5.5 


5.4 


OAO 


4948.9873 


-0.8 


4.7 


OAO 


4971.0609 


-4.8 


9.8 


BOAO 


4984.0086 


-34.5 


3.7 


OAO 


4988.9737 


-43.1 


4.1 


OAO 


4994.0597 


-29.7 


9.3 


BOAO 


5025.9988 


-9.8 


7.3 


OAO 


5026.9737 


8.0 


6.8 


OAO 



16 



Table 3. Orbital parameters of HD 100655 b 



Parameter Value 



Ki (m s -1 ) 


35.2 ± 2.3 


P (days) 


157.57 ± 0.65 


e 


0.085 ± 0.054 


u> (deg) 


132 ± 37 


T (JD) 


2453072.4 ± 15.9 


ARV* (m s" 1 ) 


-28.1 


rms (m s _1 ) 


11.2 


Reduced \[y? 


1.6 


N obs 


45 


aisinip (10~ 3 AU) 


508+ ' 034 
u.ouo_ 039 


Mm) (10- 7 Af o ) 


0071 + 0014 
u.uu/ 1_ 0014 


M 2 sim p (Mj) 


1 7+ ' 1 

i -'-0.2 


a (AU) 


7fi+ - 02 



'Offset between OAO and BOAO velocities. 
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